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Abstract 

The theories defined by Lagrangians containing second time derivative are 
considered. It is shown that if the second derivatives enter only the terms mul- 
tiplied by coupling constant one can consistently define the perturbative sector 
via Dirac procedure. The possibility of introducing standard canonical variables 
is analysed in detail. The ambiguities in quantization procedure are pointed out. 
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I Introduction 



Higher-derivative theories were introduced quite early in an attempt to regularize the 
ultraviolet divergencies of quantum field theories pQ. 

Another contex in which higher- derivative and nonlocal theories appear naturally 
is the description of the low-energy phenomena in terms of effective action which is 
nonlocal as a result of integrating out high energy degrees of freedom [2]. Moreover, 
theories with infinite degree derivatives do appear in the framework of string theory 
[3], [I] and as a modified theories of gravity [5]. 

In recent years the emergence of noncommutative field theories [6] has revived the 
discussion concerning higher-derivative theories. Apart from their string-theoretical 
origin noncommutative field theories can be viewed as an attempt to describe the 
dynamics at the scales where the very notion of space-time point lacks its meaning. 
Such theories, when modelled with the help of commutative space-time endowed with 
star product, lead at once to nonlocal Lagrangians. If noncommutativity involves 
time variable the theory becomes nonlocal in time and is plagued with unitarity and 
causality problems, at least when quantized with the help of naive Feynman rules [7]. 
There exist alternative quantization schemes which seem to cure the untarity problems 
[8]; however, they are claimed to lead to new troubles [9]. 

In view of this state of art it seems necessary to reconsider the quantization problem 
starting from first principles. First step is to put the theory in Hamiltonian form. The 
relevant framework is provided by Ostrogradski formalism [10] for higher derivative 
theories and its sophisticated version [TT] for nonlocal ones. The main problem with 
such procedures is that the resulting Hamiltonians are necessarily unbounded from 
below due to their behaviour at the infinity of phase space. This implies that the 
quantum theory, if exists, has no stable ground state. 

Still, some hope exists because, in most interesting cases, the nonlocality enters 
only throught interaction term. Then one can pose the problem of quantizing the 
perturbative sector of the theory [I], [12]. The initial value problem for perturbative 
solutions involves basic variables and their first time derivatives so the phase space 
for such solutions resembles the standard one. Moreover, for perturbation theory only 
the vicinity of the phase space is relevant and we can hope that Hamiltonian is here 
bounded from below leading to stable perturbative vacuum. 

In the present paper, inspired by Refs. [I] and [12], we study in some detail the 
Hamiltonian formalism and quantization for simple system described by the Lagrangian 
containing second time derivative in the interaction term. In Sec. II we show in full 
generality that the perturbative sector of our theory can be described with the help 
of Dirac method. There are two constraints of second kind which allow to eliminate 
perturbatively Ostrogradski momenta in favour of coordinates q and q. Dirac bracket 
{q,q}o can be then perturbatively computed to arbitrary order in coupling constant. 
It is, however, rather complicated. Therefore, in Sec. Ill we analize the possibility of 
simplifying the form of Dirac bracket and the Hamiltonian. We show that it is indeed 
possible to define perturbatively, order by order, the new variables x,x such that: (i) 
the Dirac bracket takes standard form {x,x} = 1, (ii) the Hamiltonian is the sum 
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of kinetic and potential energy. We show that there is a large freedom in defining x 
and x obeying (i), (ii); in fact, at any order of perturbative expansion for x one can 
add many terms with new, also dimensionful, constans. These constants are spurious 
in the sense that they disappear after coming back to the original dynamical variables. 
However, this might be not the case in quantum theory as we explain in Sec. VI. Sec. IV 
is devoted to the special case of homogeneous (monomial) potentials. The form of the 
transformation (q, q) — > (x, x) is studied in some detail. In particular, it is shown that 
if the degree of homogenity is odd the above transformation can be chosen such that 
the resulting Hamiltonian is parity invariant. This implies that the initial theory, when 
restricted to the perturbative sector, posses some complicated discrete symmetry. The 
form of symmetry transformation can be determined, order by order; however, we would 
like to have simpler and more straightforward explanation of its emergence. In Sec. V 
we study the simplest example of homogeneous potential of third degree, considered 
already in Refs. p[j and [12J. We find explicitly, up to fourth order, the transformation 
relating q and x as well as the Hamiltonian to this order, expressed in terms of x, x 
variables. It appears that the resulting parity invariant potential is positive term by 
term, up to fourth order. On the other hand, the initial Hamiltonian, considered to 
the same order, is not positively definite. There is no contradiction here because our 
expansions are at best asymptotic and valid at vicinity of phase space. Moreover, we 
do not know whether the property of positivity of parity invariant potential persists in 
higher orders. If this the case, the theory is perturbatively stable. 

Finally, in Sec. VI we study the quantum theory of the system described in Sec.V. To 
this end we consider (up to the second order) the transformation converting the Dirac 
bracket and the Hamiltonian into the standard form. As it has been already stressed 
such a transformation is not uniquely defined. We consider a one - parameter family 
of transformations and show that for different values of the parameter the resulting 
quantum theories are not equivalent. Specifically, the energy eigenvalues differ by an 
overall constant. Therefore, the additional parameter, spurious in the classical case, 
becomes meaningful when quantum corrections are taken into account. This shows 
that the quantization is subtle and ambiguous procedure. 



II Hamiltonian formalism for the perturbative sec- 
tor 

Let us consider the following Lagrangian 

L = \q 2 -^q 2 -gV{qA,-q), (1) 

It depends on second derivative q; however, q enters L only throught V which, in turn, 

is multiplied by the coupling constant g. 

The corresponding Euler-Lagrange equation reads 

.. 2 [dV d fdV\ d 2 fdV\\ n 
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This is a fourth-order differential equation. In order to obtain a unique solution one 
has to impose the initial conditions involving q and its first three derivatives. Corre- 
spondingly, the phase space of the system must be four dimensional. 

The canonical formalism for our system can be introduced according to the Ostro- 
gradski prescription [10J. To this end we define the canonical variables 
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dq dt\ dq 

= P2(q,<i,q) 

jq uq 

and the Hamiltonian 

H = piq 2 +P2?(gi,92,P2) - L(qi,q2,q(qi,q2,P2)) 



pi{q,q,q,<i) 



(4) 



where q(qi,q2,P2) is the solution to the last eq.Q. 

The main disadvantage of H is that pi enters it linearly so it is unbounded from 
below; the system is unstable. One can try to cure this by imposing constraints con- 
fining the system to some submanifold of phase space. A natural choice is to consider 
only perturbative solutions to eq.(T5]). Due to the fact that third and fourth derivatives 
enter only the terms multiplied by the coupling constant, the perturbative solution 
is uniquely determined by imposing the initial conditions on q and q. In particular, 
higher derivatives can be expressed in terms of q and q. In fact, one can write [12] 



q 
q- 



f{q,q) 
d_ 

dq 



i— + f(i,<i)ii)f 
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dq) 



Df 



(5) 



>) 



D n-2j 



The form of f(q, q) is determined by demanding that it is consistent with Euler- 
Lagrange equations. Let F(q,q,q,...,q^ n ') be any function. Define [12] 



[F](q,q)=F(q,q,Df,...,D n - 2 f) 

Some properties of the bracket [ • ] are discussed in Appendix. 
The consistency condition for / reads 



q+u q+g 



dV _ d fdV\ d 2 (dV\ 
dq dt\dq ) dt 2 \ dq J 



Assume now that we have found some f(q,q) obeying (J7j). 
can be now converted into constraints 

tpi=Px~ \Pi(q, q, q\ 9)](qi, qi) = 
^2=P2~ \p2(q,q,q)](qi,q2) = 



(6) 



= (7) 
The definitions of p± 2 
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Differentiating (fx, tp 2 with respect to time and using ([7j) and (jSJ) we find that there are 

no secondary constraints. 

The constraints (pi, tp 2 are second class ones: 



/ 1 9 



dV 
dq 
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_d_ 

dq 2 



dV_ d_ fdV_ 
dq dt \ dq 



(9) 



Due to the form of constraints ifi, the momenta Pi can be expressed in terms of q± and 
q 2 which parametrize the reduced phase space. Dirac bracket reads 

{A, B} D = {A, B} + {A, i Pl }{i Pl , ip % Y x {ip 2 , B} - {A, <p 2 }{ip u B} (10) 

In particular 

{qi,q 2 } D = -{tfu^y 1 (11) 
The same result is obtained by considering the symplectic form 

£1 = dpi A dq x + dp 2 A dq 2 (12) 
reduced to our submanifold. Indeed, we find 

'<9[pi] d[p 2 ] 



n 



red 



so that 



dq 2 dq x 

d\pi] d[p 2 



dq 2 A dqx 



red 



dq 2 dqx 



(13) 



(14) 



which, by eq.®, coincides with (jlip . One can also check the validity of Hamiltonian 
equations. It is convenient to come back to initial notation qx = q, q 2 = q. Simple 
computation gives 



d[H] 



f(q,q){q,q}D 



dq 

d[H] -r - x -x 
- Bf = Q{q,9} D 

where eq.(j7J) has been used. Now, first Hamiltonian equation 

q = {q, [H]}d 
gives the identity q = q while the second one 

q = {q,[H]} D 

leads to contraint equation 

q = f(q,q) 



(15) 



(16) 



(17) 
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Ill Simplifying dynamics 



The form of reduced dynamics presented above is rather complicated; in particular, 
due to the nontrivial form of basic Poisson (Dirac) bracket (lllj) the quantization poses 
nontrivial ordering problem. In order to avoid this problem one can adopt the following 
strategy [I], [12]: instead of direct quantization one first makes Darboux transformation 
which simplifies Q r ed, &>red = dx A dx. Such a transformation is not unique; in fact, it 
is defined up to a canonical transformation. The question arises whether this freedom 
can be used to simplify also the Hamiltonian or even to put it in standard form: kinetic 
plus potential energy. 

In order to analyse this problem we start with the lowest order approximation. Let 
us first note that for the Lagrangian ([!]) the zeroth-order approximation to f(q, q) reads 



fo{q,q) = -u 2 q 



(19) 



The corresponding approximation to the time-derivative operator D will be denoted 
by -Do, 



D = qTr-u q— 
dq dq 



(20) 



Finally, [ • ] denotes [ • ] given by eq.© with D replaced by D . 

Our aim is to define the transformation (q, q) — > (x, x) simplifying both Dirac bracket 
and Hamiltonian. To the first order in g one can write 



q = x + gm(x, x) 



(21) 



or 



To this order we have also 



x 



q - gm(q,q) 



(22) 



x = q - gD Q m(q, q) 



(23) 



or 



q = x + gD m(x, x) 



(24) 



where D Q on the right-hand side of eq. fl24l is given by eq. fl20l with q replaced by x. 
We start by writing the reduced symplectic form to the first order in g: 
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We are looking for m(x, x) such that the transformations ( 12TI) . ( 1241) lead to ft re d(i) — 
dxAdx. As a result of simple computation we obtain the following equation for m(x, x) 



d(D m) dm d 
dx dx dx 
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(26) 
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which we rewrite as 



d_ 

dx 



Eq.flzTJ) implies that 
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(2? 



for some function $. By virtue of (128]) $ obeys 



dx 



dV 

dx 



d[V} 



J o 



dx 



or 



D $ = [V] - V{x) 



(29) 



(30) 



V(x) being an arbitrary (up to now) function of x alone. 

In order to answer the question whether we can always find, to the first order in g, 
the transformation which puts ft r ed(i) i n Darboux form let us note that we are looking 
for a transformation which, up to a given order, is defined globally in the phase space 
(optimally, m(x,x) is some polynomial provided V is). 
Let us introduce the polar coordinates 



x = r cos 6 
x = ur sin 6 



(31) 



In terms of new coordinates eq. (l30|) reads 

_ = -uj([V] (r cos9, a;rsin9) - V(rcos6)) (32) 
The right-hand side is some periodic function of G. Therefore, one has 

|| = Y,Mr)e me + ^)e~ me ) + a (r) (33) 

n>0 

and $ is globally defined (periodic) provided ao(r) = 0. Consider the first term on the 
RHS of eq. (l32l . It is easy to see that the 9 -independent term must be a function of 
r 2 . Consider particular contribution of the form otkr 2k ; it can be cancelled by the term 
a kj2k\X 2k entering V(x). We conclude that V(x) can be chosen in such a way that no 6 

-independent term appears on the RHS of eq. (!32|) . With such a choice m(x, x), defined 
by first eq. (1281) . defines the transformation leading to standard symplectic form. Let 
us note that there is a considerable freedom in the choice of V(x). 
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In order to find the meaning of V(x) let us note that eqs. fl28|) and fl30l) imply the 
following identity 



(D% + tu 2 )m + 
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(34) 



Now, by computing x from eq. (j23|) . keeping terms up to the first order and using eq.(!33~ 
we arrive at the equation of motion for x: 



2 dV(x) 
x = —uj x — g — - — 
ox 



(35) 



Therefore, due to {x, x} = 1, the Hamiltonian computed to the first order in g, has 
the form 



H 



1 9 uj 2 x 2 
-x 2 H 

2 2 



+ gV{x) 



(36) 



Let us generalize our analysis to arbitrary order in g. To this end we write 

<x> 

q -J2g n m n (q,q) = q-M(q,q) (37) 

n=l 
oo 

9 n Dm n (q, q) = q - DM(q, q) 



x 



x 



n=l 



Let us note that the second formula does not represent an explicit expansion in powers 
of coupling constant g. This is due to the fact that D itself contains f(q, q) which is 
also given as power series in g. 

Now, assuming that Q re d takes the standard form when expressed in terms of x and 
x, we can write 



red 



dx A dx 



dDM dM dDMdM dDM dM\ , 

+ ^^--WHi) iqAiq (38 » 



dq dq dq dq 
By virtue of eqs.Q, (TT3|) and (|38|) we find that M obeys 
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(39) 



We want to solve eq. fl39|) perturbatively in g. Assume it holds up to n-th order and 
consider the n + 1 -st order. Note that the expression in the parenthesis is to be 
computed to n-th order only. Moreover, noting that M and DM are both at least 0(g) 
we conclude that the equation for n + 1 -st order contribution to M reads 



d(D m n+1 ) dm n+1 



dq 



dq 



sum of known terms 



d 2 R 



n+l 



dq 2 



(40) 



S 



where the known RHS we have rewritten for further convenience as a second derivative 
with respect to q (which is always possible). 
Eq.( )4*0|) can be written in the form 

Of dR n+1 \ dm n+1 

Again we conclude that 

m n+l = —q— (42) 

dR n+1 d<S> n+1 

D m n+1 — - = — 

oq oq 

for some 3> n+ i(q,q). Eqs. (j42l lead to the consistency condition for $ n+ i. 

D d$ w+1 <9$ n+1 = dR n+l 
dq dq dq 

or 

D $ n+1 {q, q) = R n+ i{q, q) + S n+1 {q) (44) 

One can repeat the arguments used in the case of first order approximation. Namely, 
$ n+ i is globally well-defined provided S^+i is chosen in such a way that no 9 - 
independent term (cf. eqs. (!3T|) ) appear on the RHS. This is always possible so we 
conclude that one can construct the standard canonical variables defined globally to 
arbitrary order in g. 

Let us further note that the transformation (q, q) — > (x, x) of the phase space is 
defined in such a way that the second canonical variable continues to be the time 
derivative of the first one (for a given perturbative dynamics). Therefore, the first 
Hamilton equation is an identity which, due to {x, x} = 1, leads to the standard form 
of the Hamiltonian, 

H = ix 2 + \?x 2 + V(x; g) (45) 

This can be also checked explicitly. 

We have shown that, order by order, one can reduce to the standard form the 
perturbative sector of the dynamics defined by the Lagrangian ([T|). 

IV Homogeneous potentials 

Let us now consider the special case of homogeneous monomial potentials 

V(q,q,q) = q k q l q m , m > 2; (46) 
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let us denote a = k + I + m. For dimensional reason one can write 

oo 

/M) = j> B /»(g, ?) (47) 

n=0 

where f n (q, q) are homogeneous polynomials of degree n(a — 2) + 1. 

Also, one can write the perturbative expansions for other relevant quantities. First, we 

have 

oo 

tired = (1 + Yl 9 n ^n(q, q))dq A dq (48) 

71=1 

where u> n (q, q) are homogeneous polynomials of degree n(a — 2). On the other hand, we 
have seen in the last section that there is a large freedom in the choice of the functions 
m n (q,q). Indeed, they are determined by the choice of S n (q) (cf. eq.flEJ)). There is 
only one condition restricting the admitted form of S n (q): the sum on the RHS should 
not contain the G -independent term. This is rather weak condition which allows to 
add many terms (say, any homogeneous polynomial of odd degree) containing new 
(also dimensionful) parameters. However, one can show that it is always possible to 
choose the " minimal" S' n s in the sense that the only constants entering them are g and 
oo. Assuming this is the case up the order n we conclude that 9 is homogeneous 
polynomial of degree (n + l)(a — 2) depending only on one constant uj. Therefore, 
R n+ i can be also chosen as homogeneous polynomial of degree (n + I) (a — 2) + 2 
containing only one dimensionful constant u. As a result, the -independent term 
in Rn+i must be of the form r ( n + 1 )( a ~ 2 )+ 2 times a dimensionless constant. Then we 
can choose S n+ \{q) as proportional to g( n+1 K a ~ 2 )+ 2 and $ n +i obeying eq.(144"l) can be 
taken as homogeneous polynomial of the same degree depending only on uj. So, by 
first eq. pSl) m n+1 is homogeneous of degree n(a — 2) + 1. This concludes the inductive 
proof. 

With the minimal choice of the transformation ( l37j) one can easily write out the 
general form of the potential V(x; g) entering the Hamiltonian (T4"5l) ; it reads 

oo 

V(x; g) = J2 v n g n uj {l+2m - 2)n+2 x ia ~ 2)n+2 (49) 

n=l 

Let us now consider the particular case of odd a. Notice that R n is of degree n(a — 2) + 2 
which is odd for n odd. Therefore, R n is then homogeneous polynomial of odd degree 
so it does not contain B -independent term. So S n can be chosen as a n q n<ya ^ +2 with 
a n arbitrary (in particular, one can take a n — ). It is not difficult to see that a n can 
be chosen perturbatively order by order so that the odd terms in the expansion fT4"9"l) 
vanish. Indeed, let 

S n = aq n{a ' 2]+2 (50) 

Once S n is selected, one can define, via eqs. (1371) . fl4"2"l) and (JHJ, the variables x a ,x a 
to n-th order. It is easy to see that the relation between x ,xo (corresponding to the 
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choice a — ) and x n , x„, to the same order, reads 



with 



x = x a + g n A m a (x a , x a ) (51) 
x = x a + g n D A m a (x a , x a ) 



9Af n „ A -9A$ n 

Am ft = — — — , D Am a = 

ox ox 

D A $ n = ai n(fl - 2)+2 (52) 



Therefore, adding the term ([50]) amounts to the following change of the Hamiltonian 

tt • 2 i ^ 2 2 i t> / \ 1-2 , 1 22 , 

2 X ° 2^ X ° n ^ ' ^ ~ 2 Xa + 2 UXa 



+V n (x a , g) + g n x a D Am a + g n u 2 x a Am, 

d 

dx sy dx 



\x 2 a + \uj 2 x 2 a + V n (x a , g) - g n ( x a -^- u 2 x a -^- ) A $„ 



~^ + + V n {x a ,g)\ - a« (a " 2)+2 (53) 

Adjusting properly a one can cancel, order by order, all odd terms in V(x, g). 

Concluding, we find that for odd monomial V(q, q, q) one can reduce, order by 
order, the perturbative potential V(x; g) to the form 

oo 

V(x- g) = J2 v 2k gV +2m - 2 ^ k+2 x 2 ^ k+2 (54) 

k=l 

Note that in this case the perturbative sector exhibits some discrete nonlinear symme- 
try. In fact, the resulting standard Hamiltonian is parity invariant: X ^ X . X ^ X 
is a symmetry. Then, expressed back in original variables, the parity transformation 
produces nonlinear symmetry defined order by order in coupling constant g. 



V The simple example 

Let us consider a simple model studied already in Refs. [4], [12] : 

r *2 1 2 2 -2 / cc \ 

L = -q --wg -gqq (55) 

It belongs to the class of models studied in the last section. Eq. ([55l) leads to the 
following equation of motion 

q + u 2 q + g{3q 2 + 4q q +2qq (IV) ) = (56) 
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The canonical variables read 



Qi = Q, Q2 = q 

P! = q + 2g{qq + q q) (57) 
P 2 = -2gqq 

It is also straightforward to write out the Hamiltonian 

P 2 1 1 

H = P x q 2 - ql + -u?q\ (58) 

y Agq x 2 H2 2 yi v ; 

In order to perform the reduction to the perturbative sector we impose the constraint 

q = f(q,q) (59) 
Then, by virtue of eq. f[5"6l) . f(q, q) obeys 

f + Uj2q + g { 3f2 + M {% + f %)) + m 

+2? ? j2+ IJTir + fjT + ^TTF 7 + / F +^7F^ =0 
V oq z oq oq oq Oqoq \oq J Oq z J 



This equation, although quite complicated, can be solved perturbatively order by order 
in g. For example, to the third order in g one finds 

/ = -u 2 q - g(hu A q 2 - Auj 2 q 2 ) + # 2 (-76aA 3 + UOco^qq 2 ) + (61) 
-# 3 (1959u;y - 6800cAV + 736u; 4 g 4 ) 

The constraints (jSJ) take the form 

Pi- (7-20 (qf + qM + qM) « 

F 2 + 2gqf ~ (62) 

while fi re d is given by 

fired = (63) 

<9f 9f d 2 f fdf\ 2 d 2 f\ 

l + Agf + Agqjf- + 2gq-^ + 2gqq—j- + 2gq[-±) + 2gqf—^ ) dq A dq 
oq oq oqoq \oq J oq z I 

Finally, the reduced Hamiltonian reads 

[H] = \q 2 + \u 2 q 2 + g (~qf + 2q 2 f + 2qq 2 ^ + 2qqf^) (64) 
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Now, one can try to find perturbatively the "normal" coordinates x,x. Following the 
method outlined in previous sections we found that, to the fourth order, 




(65) 



and 




30136 



45 



gW°x 6 + 0(<? 6 ) 



(66) 



We see that our perturbative Hamiltonian, when put in normal form, becomes posi- 
tively defined, at least up to fourth order in g. We don't know whether this property 
persists in higher orders. Let us note that our reduced Hamiltonian (151)1 is not positive. 
For example, to the first order in g one finds from (I6T!) and (1641) 



which is negative for large q, q. 

On the other hand, to the same order [H] , when expressed in terms of new coordinates, 
is simply the energy of harmonic oscillator. We conclude that, at best, we can expect 
that our series defining new coordinates are asymptotic (note that [H], as given by 
eq. flBTl) . becomes negative for q, q of order ^ ). 

VI Quantum theory 

Our ultimate goal is to quantize the higher derivative dynamical system. The main 
disadvantage of the Hamiltonian formalism introduced by Ostrogradski is that some 
momenta enter the Hamiltonian linearly. Therefore, it is unbounded from below. Con- 
trary to the case where the Hamiltonian is unbounded in small regions of phase space, 
this kind of unboundness cannot be cured with the help of uncertainty principle. As a 
result, no stable ground state can exist. 

However, one can ask whether it is possible to quantize consistently the higher- 
derivative theory in the perturbative sector. The first trouble is related here with the 
complicated form of reduced symplectic structure. It is by far not sure whether one 
can find the proper ordering procedure which allows to convert complicated Poisson 
brackets into commutators obeying Jacobi identity. 

The simplest way to define the perturbative quantum theory seems to be the fol- 
lowing. First, we construct on the classical level the transformation in reduced phase 
space leading to the standard from of the Poisson bracket and the Hamiltonian. Then 
the quantization can be performed in a straightforward way. Moreover, if the classical 




(67) 
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Hamiltonian appears to be bounded from below, the quantum theory possess pertur- 
batively stable ground state. Once the theory is quantized in "standard" coordinates 
one defines the quantum counterparts of initial variables by inverting (perturbatively) 
the classical map and choosing a definite ordering (for example, the Weyl one). 

The main problem here is that such a procedure is by far not unique. In fact, we 
have seen in previous section that there is a large freedom in defining the classical trans- 
formation to standard coordinates. One can hardly believe that the quantum theories 
resulting from different choices of such transformations are equivalent. Moreover, in 
the process of defining the perturbative transformation from q - to x -variables one can 
introduce new (also dimensionful) constants. On the classical level they are spurious 
and disappear after coming back to original dynamical variables. This may be not the 
case after quantization has been performed and the additional parameters may appear 
to be relevant. 

In order to illustrate this phenomenon let us go back to our simple model. Consider 
the transformation 

x = q + g{(3u 2 q 2 + (2(3 + 4)g 2 ) + 

g\(-2(3 - ^-)u 4 q 3 + (-32 - 2(3 2 - 2A(3)u 2 qq 2 ) 

x = q- 2g((3 + 4)u 2 qq + g 2 ((4(3 2 + 22/9 - 2Q)u A qq 2 - 2((3 2 + 4{3)uj 2 q 3 (68) 

depending on one real parameter (3. In terms of new variables the Hamiltonian takes 
the form 

[H] = h 2 + ±wV - g{(3 + 1) W V + 5g 2 Q/3 2 + (3 + ^ ^ 6 x 4 + 0(g 3 ) (69) 
For (3 = — 1 we obtain the parity invariant form. 

Let us now compute the energies to the second order in g. Standard perturbation 
theory gives 

E n = f^(n+ 1 -^ + 2 -^g 2 h 2 u\n 2 + (n + l) 2 ) + \g 2 h 2 uj\(3 + l) 2 (70) 

We see that the energy eigenvalues depend on (3, although it is only an overall shift. 
It is interesting to note that the energies take minimal values in the parity - invariant 
case. 

The ambiguity considered above is rather mild. We could add other terms, much 
more complicated and containing new dimensionful constants. Let us remind that 
the only condition imposed, order by order, on new S n (cf. eq. (l4"4"|) ) is that the RHS 
contain no O -independent terms. Keeping this in mind one can easily understand that 
the resulting form of standard Hamiltonian can vary considerably depending on the 
particular transformation chosen. This may have strong impact on the form of energy 
spectrum. The resulting quantum theories become nonequivalent. This effect can be 
ultimately ascribed to the ordering problem. 
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VII Appendix 



The main property of the symbol [ ■ ] introduced in Sec. II is expressed by the equation 



d[F] 
dt 



D[F] 



dF 
~dt 



(71) 



To see this let us note that 112 



d\F\ 
dt 



Due to [q — f] = we find 



dF' 
~dt 



\d[FY 




'dF' 


. dt . 




. dt . 



Also D[F] = [D[F]] and D[F] = [D[F]} 
Iterating ( J7TT) one obtains 



d[F] 
dt 



[ dt ] 



'd 2 F~ 




' d 


'dF' 






' d 


\d[FY 






'd 2 [F}' 


dt 2 




dt 


dt 






di 


dt 






dt 2 



and 



D 2 [F] = D[D[F}} 



dt 



(D[F)) 



d 2 [F] 
dt 2 



Therefore 



'd n [F\ 
dt n 



D n [F] 



'd n F' 
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